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§ '. Abstract 

A theory of bipolaron states in quantum wires with a parabolic potential 
\J3 ■ well is developed applying the Feynman variational principle. The basic pa- 

rameters of the bipolaron ground state (the binding energy, the number of 
phonons in the bipolaron cloud, the effective mass, and the bipolaron radius) 
are studied as a function of sizes of the potential well. Two cases are con- 
sidered in detail: a cylindrical quantum wire and a planar quantum wire. 
Analytical expressions for the bipolaron parameters are obtained at large and 
small sizes of the quantum well. It is shown that at R S> 1 [where R means 
the radius (halfwidth) of a cylindrical (planar) quantum wire, expressed in 
Feynman units] , the influence of confinement on the bipolaron binding energy 
is described by the function ~ 1/R 2 for both cases, while at small sizes this 
O ■ influence is different in each case. In quantum wires, the bipolaron binding 

!> ■ energy W (R) increases logarithmically with decreasing radius. The shapes 

and the sizes of a nanostructure, which are favorable for observation of stable 
bipolaron states, are determined. 

I. INTRODUCTION 

Landau's idea [Q] of the auto-localized state of a charge carrier (polaron) in a homoge- 
neous polar medium got a further development by Pekar who first studied a problem of 
a stable complex of two charge carriers of the same sign (bipolaron). The bipolaron binding 
energy was first calculated in Ref. |J. The bipolaron problem was widely discussed, see e. g. 
Refs. @-§|. A detailed outline of this subject is presented in the recent review |T0[| . 
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Dimensionless constants of the Coulomb interaction U and of the electron-phonon inter- 
action a are related to each other by the equation |LQ| : 



U = ^, (1) 
1 — 77 

where rj = £00/^0 ( £ o an d e oo are static and optical dielectric constants, respectively). Due to 
the fact that e > e^, the relation U > \[2a follows. When the distance I between electrons 
is large or small compared with the characteristic polaron radius R p (see Ref. ||11||), the 
phonon-mediated attraction between electrons occurs to be weaker than the repulsion. At 
large distances I 3> R p , both interaction potentials have similar spatial dependences but the 
Coulomb repulsion is stronger than the phonon-mediated attraction. In the opposite case, 
I <C Rp, the Coulomb potential diverges at the zero distance, while the phonon-mediated 
attraction is always finite. Nevertheless, when two electrons move in such a way that the 
average distance between them is of the same order as the polaron radius, the bipolaron can 
be stable at a ^> 1 and When two electrons are together confined to a potential well, 

one can expect that the conditions of the bipolaron stability may be improved for relevant 
sizes of the well. 

Two new circumstances have stimulated the bipolaron theory: the progress in the fab- 
rication technology of mesoscopic nanostructures such as quasi-2D (quantum wells and su- 
perlattices), quasi-lD (quantum wires), quasi-OD (quantum dots), and the advancement of 
the hypothesis that bipolaron excitations might play a role in processes occurring in the 
high-temperature superconductors. The present research has been motivated also by the 



recent advances in creation of nanocrystals with a strong ionic coupling [12 . 

The basic bipolaron parameters are recalled in what follows. The bipolaron stability 
region is determined by the inequality W > for the bipolaron binding energy 

W = 2E P — E hip . (2) 

Here E p and E^ p are the free polaron and bipolaron ground state energies, respectively. 
From the equation 

W(a, V ,TZ) = 0, (3) 

where TZ denotes the set of parameters determining the shape and the size of the confinement 
domain, the functions a c (77, TZ) and r\ c (a, TZ) describing the boundaries of the bipolaron 
stability region are found, for fixed 77 and a, respectively. According to different theoretical 
treatments PHITiH the bipolaron binding energy is an increasing function of a and a decreasing 



function of rj. It will be shown that the function rj c (a,TZ) starts from r] c = at a = 
awn (TZ) 7^ 0, grows with increasing a and tends to the upper limit ^ max at a — >• 00. 
The bipolaron stability region is then determined by the inequalities a > ct m i n (TZ) and 
< rj < rj c (a,TZ). 

Let us adduce typical values of the parameters a m in,3D an d ^max,3D of the bulk (3D) 
bipolaron: a mirii 3D = 6.8 and r/ ma x,3D = 0.14 were found by Verbist, Peeters and Devreese 
]T3] , P^ | and by Verbist, Smondyrev, Peeters and Devreese Adamowski [^] obtained 

tt m in,3D = 7.3 and ^ maXj 3D = 0.14. The bipolaron theory developed for pure 2D fl6UTT|l and 



ID [17| models shows that the bipolaron stability region broadens when the dimensionality 



2 



is reduced. For these systems, the following parameters were obtained: a m i nj 2D = 2.9, 
^max,2D = 0.158 (Ref. fL6|); a mirijlD = 0.9, r? maXi i D = 0.764 (Ref. fl8|l ). Bipolaron states 
were investigated in a quantum well [|19],|20| and in a quantum wire EH] as a function of the 
characteristic size of the system. The polaron theory for a quantum dot is developed in 



Refs. 22-25 



The goal of the present investigation is to determine the bipolaron stability region and 
to study the basic parameters characterizing the bipolaron ground state as a function of 
confinement. Two different types of confinement are considered and compared to each 
other: (i) a cylindrical quantum wire of the radius R, where continuous transitions from 
3D to ID are realized with decreasing R; (ii) a planar quantum wire of the width L, where 
a transition from 2D to ID is realized with decreasing L. A unique approach, namely the 



Feynman variational method |26|,[^]], is used throughout the paper for both systems under 
analysis. 

The paper is organized as follows. In Section II, general formulae for parameters of a 
bipolaron in quantum wires are deduced. In Section III, particular cases of cylindrical and 
planar quantum wires are considered. The basic parameters of the bipolaron ground state 
are obtained. Limiting cases of strong and weak confinement are studied in detail. The 
obtained numerical and analytical results are discussed in Section IV. Section V contains 
conclusions about the influence of confinement on the bipolaron binding energy in quantum 
wires. 



II. GENERAL THEORY 

We analyze the bipolaron problem taking into account both the electron-phonon inter- 
action and the Coulomb repulsion between two electrons confined to a quantum wire. The 
Lagrange function of the system is 

miX A in e 1 



E E E «< 

4=1 n=l,2 Z n=l,2 - ^ 



+ 9 E («£ - "o w l) - E E 7k (r„) w k , (4) 

Z k 71=1,2 k 

where r n (xi n , X2> n , x^ n ) is the radius vector of the n-th electron (n = 1,2); rrii is the ii-th 
component (i = 1, 2, 3) of the diagonal band mass tensor, U (r) is the potential energy 
of an electron in the quantum wire, are the normal coordinates of longitudinal optical 
(LO) phonon modes. Here, the parameter D determines the dimensionality of the electron 
subsystem: D = 3 and 2 for cylindrical and planar quantum wires, respectively. Amplitudes 
of the electron-phonon interaction are taken in the Frohlich form: 

where fh = {mim2m^) 1 ^ i , V is the volume of the system, and the Frohlich constant 

a= _fL(j_-iu^y /2 (6) 
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characterizes the strength of the coupling between an electron and bulk polar LO phonons 
with the long- wavelength frequency u Q . In this paper, the 3D phonon approximation is used, 
according to which the interaction of an electron with both bulk-like and interface phonons 
is replaced by that with 3D phonons. This often used approach is adequate because any 
integral polaron or bipolaron effect, resulting from a summation over all phonon modes, 
appears to be only weakly dependent on the details of the phonon spectrum. It should 
be also mentioned that the system under consideration simulates realistic structures with 
relatively smooth interface barriers, where interface-like phonon modes can appear, which 
are smoothly distributed in space rather than localized near a sharp boundary, as is the case 
for interface modes. 

In order to study the bipolaron problem at arbitrary values of a, the Feynman variational 
approach [[T^J is the most appropriate method. The trial Lagrange function is written as 



1 

Ltr = — [ m *^?,n MiXi n — hi (Xi t n — Xi jTl ) — h\ [Xi >n — Xi^, 

^ i=\ n=l,2 

+ (xa - x t2 f - ]T W (r ft ) , (7) 



f 

i=l n=l,2 



where Xi n are coordinates of the n-th "fictitious" particle (n = 1,2). This model imitates 
the interaction of electrons with phonons and between each other by elastic bonds as shown 
in Fig. 1. The masses Mj and the force constants hi, k' { , Ki play the role of variational 
parameters. For n — 1, n takes the value 2, and for n = 2, n is equal to 1. The potential 
well U (r) from Eq. (^) is simulated here by a parabolic function: 

W(r) = i (8) 

Z 8=1 

The index q characterizes the dimensionality of confinement and is determined as follows: 
for a planar quantum wire q = 1 (Qi ^ and Q 2 = 0) and for a cylindrical quantum wire 
q = 2 (fij ^ 0, tt 2 ^ 0, and fi 3 = 0). 

The basis of the Feynman variational method is the Jensen-Feynman inequality Wt 



(exp (S - Str)) Str > exp (S - S tr ) Str , (9) 
where the angular brackets denote averaging over electron paths: 

frs _ Tr/DrG'[r]exp(g tr ) 

{ )s *- Tr/Drexp(,S ir .) ' 1 Uj 

Here S and S tr are the electron action functionals obtained after integration over phonon 
variables and over coordinates of "fictitious" particles, respectively. At low temperatures, 
the variational bipolaron energy is calculated using the expression 

E Up = E tr - Jim ^'f^ 5 " , (11) 

where E tr is the ground state energy of the trial system with the Lagrangian ([7|), f3 — 1/kgT 
is the inverse temperature. 
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The trial Lagrange function (0) consists of D independent parts: L tr = Li. Each part 

Li is a function of four variables xn, 2^2, -Xji, X i2 . Let us introduce unified denotations 
for coordinates of electrons and of "fictitious" particles: xn = xn, Xi2 = Xi2, xa = Xn, 
Xi4 = Xi2- It follows from the form of the trial Lagrangian ([?]) with Eq. (|8|) that the groups 
of variables Xjj with different indices i are dynamically independent from each other. They 
are related to normal variables £jj by the unitary transformation: 



i'=i 



1,...,D 



(12) 



with 4x4 matrices (j, j' = 1, . . . ,4). From the equations of motion for the group of 

coordinates Xij (j = 1, . . . ,4) with a fixed i, the following eigenfrequencies are obtained: 



LO 7 



LO, 




mi J 



Mi\ n Ki 

in; / m,- 



\ 



l-^|.?-i>: 



i 2 



M 



+ 4— v. 



mi 



3 = 1,2, 



\ 



1 M A 2 02^0^ 

1 k? - n? + 2— 

rrij J mi 



2 _|_ 4.^ ^0' 



J =3,4, (13) 

where vf = (ki + fcQ/Mj. Matrix elements of the unitary transformation (|i~2"D are 



2 _ Ul-V? d 2 
u i,H ~~ o I, .2 , ,2 V «,12 



LO 



i'2 



2 (^ - c4) ' 



2 (wfi - w i2j 



r/ 2 

u i,13 



i3 



2(u4 



r/ 2 

u i,14 



2(^3 



4)' 



»,2j' 



Sj'di t 



if, 



d 



ki -\- Sj'k; 



i,3j' 



to, 



1. 1' 



di 



d 



Sj>d ii3j >; 



(14) 



l(j 



1,2); Sj 



-1 (J = 3,4). 



Note that the elastic repulsion imitating the Coulomb interaction gives a contribution to 
the eigenfrequencies with j = 3 and 4 through the force constants Ki. It is easy to see from 



Eq. (13) that under the conditions of a strong confinement along the z-th coordinate axis 
the eigenfrequencies with j = 1 and 3 corresponding to the motion of the bipolaron along 
this whole are determined mainly by the parameter fij. 

The action functionals S and S tr in Eqs. (H) to ([0]) contain the potential energies U and 
W, respectively. Though the shape of a real potential U may differ from that of the model 
quadratic potential (||), the averaged difference (U — VV) 5tr can be omitted as far as it is 
small when compared to the rest of (S — S tr ) s //3- 

The averaging procedure in Eq. (|TT| ) is carried out by the path integration and leads to 
the following form of the variational bipolaron energy: 



trip 



J2B l + C + P. 



(15) 



i=l 
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Here the terms Bi include the averaged kinetic energies of two electrons and of two "ficti- 
tious" particles as well as the averaged potential energy of the elastic interaction of these 
four particles: 



B; 



E 

3=1 



UJij 1 - 



(jJiA — ^ 



UJ 



d i,lj 
ij / 



D. 



(16) 



In Eqs. (|I5D ) (|T6| ) and further on, the Feynman units |^S] are used: Tiujq for energies; ujq for 
frequencies; and (h/fhuJo) 1 ^ 2 for lengths. 

The averaged potential energy of the Coulomb electron repulsion is 

" r^(0), 



C (1-77) 7T 2 ' 

and the averaged energy of the electron-phonon interaction is 



00 

P = -4E / dre- T JC n (r) 

71 n=l,2„ 



where 



00 



E *? An (r) 



8=1 



(17) 



(18) 



(19) 



i=l 



The functions A, 



An (r) = — 

mi 



r. 



are determined as follows: 



E 



r7 2 



.7=1,2 ^'j 



u ij T ^j _|_ ^ 



r/ 2 



j=3,4 



n —ujijT 



1 \« 

■1) e 



n = l,2. (20) 



In order to find the bipolaron energy, it is necessary to minimize the function E^ ip given by 
Eq. ( |15D over twelve independent variational parameters uu, 



131 



1 ... 3, j = 1 ... 4, which 
are used instead of the masses and the force constants ki, k' t , Ki. 

From Eq. flllf ) for the eigenfrequencies, the expression for components of the diagonal 
tensor of relative bipolaron effective mass is deduced straightforwardly: 



rrii 



+ 1 



mi 



2(w? 1 +o;2,-n 



(21) 



where the values of parameters u>n, uj i2 , ^% are taken which provide the bipolaron energy. 
The number of phonons in the bipolaron cloud is determined by the general expression 



of Ref. ITT 



dS 



ph 



,d(huj )/ St 

which gives in the case under consideration: 



(22) 



ph = 4 E / IC n (r)e- T rdr. (23) 

71 n=l,2 { 

\) are performed using 



Calculations of the average number of phonons according to Eq 
the results of minimization of the bipolaron energy. 
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III. BIPOLARON IN CYLINDRICAL AND PLANAR QUANTUM WIRES 



A. Variational problem 



Here we write down the variational bipolaron energies in the cylindrical and planar quan- 
tum wires (see Fig. 2). Hereafter, the following denotation for the confinement parameter is 
used: Qi = Q±, i — 1, q. The electron mass is taken to be isotropic, i. e. mi = rri2 = = m. 
From Eqs. ([15]) to ([18|) we obtain the variational bipolaron energy 



E 



bip 



B± + B\\ + C + P, 



where, in accordance with Eq. fllED , 



Z 3=1 



(24) 



(25) 



B ] 



E U ±3 ( 1 b ~ 2t 'J 



i=i 



(26) 



Here the frequencies of the motion along the z-axis (called below the longitudinal motion) 
are 



u m = cjp u D2 = 0, ujbj = j = 3,4, 



(27) 



VD = v\\ 

and those of the motion in the xy-plane (the transverse motion) are 



1,...,4, v i = v± 



(28) 



with i — 1, 2 for q — 2 and i — 1 for g = 1. In the case under consideration, the integrations 
containing the function K n (r) are performed analytically. The calculation of the integrals 
in Eqs. ( |T7D and flT8"D yields the averaged potential energy of the Coulomb repulsion between 
electrons 



C 



V2U ( A ±2 (0Y 



and the averaged energy of the electron-phonon interaction 



(29) 



P 



2a 



E / drer 

1 O " 



n=l,2 



(30) 



where 



% to = E — 

j=3,4 



1 + (-l) n e" 



+ 



a D,ll 
^||1 



1T +4 u r; 
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A i 



rl 2 

3=1,2 U ±j 



-U ±j T 



d 2 

y] j 



1 + (-l)V 



n 



1,2; 



tanh 1 -y/x 
dip 



J- 

(l— a: sin 2 </,) 



T72; 



3 
Q 



(31) 



The minimization of the variational bipolaron energy Eu P determined by Eqs. fl2"4]) to (|3TD 
is carried out with respect to eight variational parameters v\\, u± 2 , w±j (j = 1,3,4). 
Setting k[ and Ki equal to zero in these formulas, the twice value of the polaron energy |25| 
is obtained from Eq. (^4j). The binding energy is then found according to Eq. (Q). Results 

— 1/2 

of the calculation of W as a function of the quantum wire radius R = Q ± are presented 
in Fig. 3 for different values of a. Then the functions a min (R) = a c (r] = 0, R) (Fig. 4) and 
T) c (R, a) are obtained. The latter function is used for calculation of the critical value of the 
Coulomb repulsion constant U c (a) in order to describe the bipolaron stability region shown 
in Fig. 5. From Eq. fl2"T|), taking into account Eq. (E7|), the relative bipolaron effective mass 
of the longitudinal motion is derived as 



(™bip). 



1*4 



ni 



(32) 



Plots of the relative bipolaron mass as a function of R are shown in Fig. 6. A detailed 
discussion of the results will be given in Section IV. 

B. Weak size quantization 



For a weak size quantization, the eigenfrequencies of the transverse motion can be rep- 
resented as expansion series in fl±. In these series, we take into account only two first 
terms: 



2 2 



, W ±2 



n 



2 V \\ 



'OJ 



0(tt 



111 



2 2 



13 



9 9 

"^3-^114 



2 2 



, ^_L4 



9 9 



As a consequence, Eq. fl26|) takes on the form 



B± =q 



Bn + n 



+ o(n] 



"4^11! Wjf! 

and the parameters in the right-hand side of Eq. (|3~0D satisfy the relations 



1 

_( 

4' 



..4 : Jr) = ..An,, ( r) - ^±r 2 ^ + 0(^1), n = 1, 2. 

^Ili 



(33) 



(34) 



(35) 



S 



Substituting the expression (|35|) in Eqs. (|29|). 
repulsion 

V2U f lq 



c 



V^II12(0) 



0), we obtain the energy of the Coulomb 



(36) 



+ o(n 2 ± ), 



and the energy of the electron-phonon interaction 



P 



2a 



'JXJ 

— 1 o 



n=l,2 



flq 



\JAUt) 



+ 



.2 V \\ 



flq 



where 



lq 



1, ? : 



2; 
1. 



2g 



9 
Q 



T) 

2; 
1. 



i), 



(37) 



The bipolaron energy in this limiting case can be represented as 

E _ pO , a p , n/n2 



>bip 



where i?° ip is the bipolaron energy in three or two dimensions for D = 3 or D = 2, respec- 
tively. The confinement- induced shift of the bipolaron energy Ai^ip is given by 



AE, 



hip 



CO 



111 



Q OO 

II 1 




It is worth mentioning, that in the strong coupling regime, the integrals in Eq. (|38 



(38) 



arc 



calculated analytically, and the minimization of this variational function with respect to the 
frequencies is performed explicitly. For this purpose we use the results of Ref. [|8[], where 
the following analytical expressions for frequencies are obtained: to\u = a 2 u)i (for % = 1,3), 
oju = l,v\\ = 1, where 



128 HW 



u 3 



128 
( 

9tt 



i-C 2 (u) 



(39) 



C(U) 



u 

16a 



D 
D 



3; 
2. 



Using these frequencies, we find the confinement- induced shift of the bipolaron energy to be 

A£bi P = q^-. (40) 
2a z u>3 



This result differs qualitatively from that deduced in Ref. |29| for a bipolaron in a weak 
magnetic field, where the cyclotron frequency uj c plays the role of Q. Namely, as distinct 
from Eq. in the equation from Ref. [[29|] for the first-order correction to the bipolaron 
energy a 2 stands instead of a 4 . 

It is important to note, that this positive correction to the bipolaron energy due to the 
confinement is less than the twice value of the respective correction to the polaron energy 
] . The confinement- induced variation of the bipolaron binding energy obeys the inequality 

9tt 1 



> 0. 



.26 D £3- 

Thus, the enhancement of the bipolaron binding takes place due to the confinement. 



(41) 
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C. Strong size quantization 



In the limiting case of a strong size quantization, the terms of the order of Q 2 ^ play a 
determining role in Eq. fll5l). For the frequencies of the transverse motion, the expansion in 
inverse powers of Q 2 gives: 



ul 1 = nl + —vl + o(ni 2 ), u 2 ±2 = vi + o(n] 



m 

^i 3 = ^i + ^i-2^ + o(fi: 2 ), ui 4 = vi + o(Qj_ 2 ). (42) 

Consequently, the bipolaron ground state energy is described by the expression 

E hip = qtt ± + 2]| + = In ( ^O ± A,| 2 (0) J 
/2ttA|| 2 (0) W J 



B i^/^CT h (^ w )- (43) 



The first term in the right-hand side of Eq. (f43|) is the energy of two electrons in the parabolic 
potential. The last three terms in the right-hand side of Eq. flilf) are due to the electron- 
phonon and Coulomb interactions. In the strong coupling regime, = a 2 Ui (i = 1,3); 
and v\\ are proportional to a with coefficients which are functions of fl±. Then omitting 
the terms of the order of a in the last three terms of the variational bipolaron energy ( [43]) 
we obtain 

In 



q 2 ot 2 uj\ 




(44) 

Note, that the second term in the right-hand side of Eq. ( f44"|) is proportional to a 2 , as is 
expected in the strong coupling regime. When replacing Q± —>■ uj c at q = 2, this polaronic 
term coincides with that of the bipolaron variational energy in a strong magnetic field from 



Ref. ||29|| . The dependence of the bipolaron binding energy on the cylindrical confinement 
provides a possibility for a controllable enhancement of the bipolaron binding by decreasing 
the radius of a quantum wire. 

Setting U = and Co = uj\ = uj 3 in Eq. (f4*4|), one obtains the twice variational polaron 
energy 2E p with the variational parameter Co (see Ref. ||25|| ). The polaron energy E p depends 
on the confinement parameter similarly to E^ ip . The binding energy W (which is not written 
explicitly to save space) increases logarithmically with increasing Q±. 



IV. DISCUSSION OF NUMERICAL RESULTS 

Beyond the framework of the limiting cases which allow an analytical treatment as dis- 
cussed above, the bipolaron stability is studied using the following computational procedure. 
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First, we evaluate the bipolaron energy E^ v and the model bipolaron effective mass defined 
as m bip = 2(M|j + m). Second, the functions a min (R) and r\ c (a, R) are found from Eq. (||). 
The region of the Frohlich coupling constant ranging from 2 to 4 is chosen for the numerical 
work in order to include the values of a corresponding to some specific substances with 



small rj (for example, Ti0 2 : a = 2.03, rj = 0.035 |3; T1C1: a = 2.56, rj = 0.133 BaO: 
a = 3.23, rj = 0.118 Q; LiBr: a = 4.15, rj = 0.24 0). 



Fig. 3 illustrates the size dependence of the bipolaron binding energy in cylindrical and 
planar quantum wires, for a values corresponding to the above-mentioned substances and for 
rj = 0. As seen from Fig. 3, the bipolaron binding energy monotonously rises with increasing 
transverse confinement [cf. Eqs. ([24]) to ([30])]. 

In Fig. 4, the minimal value a m \ n is represented as a function of R and L for cylindrical 
and planar quantum wires, correspondingly. In the region of large R and L, the minimal 
values a m in for cylindrical and planar quantum wires tend to the three-dimensional and two- 
dimensional limits a m in,3D and a mint2 D, respectively. When R and L decrease from 1.0 to 0.1, 
a rapid diminution of a m i n is seen. Note that at small values of R, L (which are, however, still 
compatible with the continuum description), the bipolaron stability region extends to small 
values of a. Note that the bipolaron parameters for quantum wires obtained in the formal 
limiting cases R,L-^0 differ substantially from those derived for the purely one-dimensional 
model |18[ (with ID-electrons and lD-phonons), which gives a minj iD = 0.9. 

In Fig. 5, the ratio of the critical value of the Coulomb repulsion constant to the Frohlich 
coupling constant a, 

= ^ (45) 
a 1 — T] c (a) 

is plotted as a function of a for various radii (ranging from 0.01 to 20.0) of the cylindrical 
quantum wire. Since the parameter r] c is non-negative, U c (a) /a cannot be less than the 
value \/2 (shown by the line A). When increasing a, the right-hand side of Eq. (|45|) tends to 
the three-dimensional limit y/2/ (1 — ?7 c ,3d), marked by the line B. The physical sense of this 
trend consists in the following: when increasing the electron-phonon coupling, the electron 
confinement to the parabolic potential (||]) is gradually replaced by the confinement to the 
polaronic potential well. The regions of bipolaron stability can exist only between the lines 
A and B. The domain between a curve U c (a) /a and the line A is the bipolaron stability 
region for a specific radius of the cylindrical quantum wire. This figure illustrates clearly 
an enlargement of the bipolaron stability region with decreasing the radius of the quantum 
wire. An analogous dependence of the bipolaron stability region on the width takes place 
for the planar quantum wire. 

Bipolaron effective masses are represented in Fig. 6 as a function of the dimensionless sizes 
R and L of the cylindrical and planar quantum wires. The size dependence of the bipolaron 
effective mass is qualitatively similar to that of the ground state energy but appears to be 
substantially more pronounced. At small radii (0.1 < R < 0.2), the bipolaron mass strongly 
increases with decreasing R. 
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V. CONCLUSIONS 



The conclusion of our analysis is that the confinement leads to an enlargement of the 
bipolaron stability region in cylindrical and planar quantum wires as compared to the cor- 
responding regions of infinite three-dimensional and two-dimensional systems, respectively. 
For R ~ 1 or L ~ 1, the critical values a c required for bipolaron stability are close to those 
for Ti02, T1C1, BaO and LiBr. For example, according to Fig. 4b, the bipolaron stability 
region sets in at the width L of the planar quantum wire of about 8 nm for parameters of 
T1C1. In this view, manifestations of the bipolaron phenomena might be observed in the 
technically achievable planar quantum wire structures. 

The performed analytical and numerical analysis of the influence of confinement on 
the bipolaron binding energy has shown that stable bipolaron states are possible even for 
intermediate values of a (a ~ 2) and for not too small values of i] (r/ ~ 0.1) in nanostructures 
whose sizes are of the same order as the polaron radius R p . Among the considered systems, 
the most favorable conditions for the bipolaron stability take place in planar quantum wires, 
where the binding energy W increases monotonously (logarithmically) with strengthening 
confinement. Nanostructures, whose sizes satisfy the conditions of the bipolaron stability, 
seem to be achievable for the modern technology. 
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FIGURE CAPTIONS 



Fig. 1. A scheme of the trial system which contains two electrons connected with two 
"fictitious" particles through the elastic attraction and models the Coulomb interaction by 
the elastic repulsion. 

Fig. 2. A scheme of cylindrical (a) and planar (b) quantum wires. 

Fig. 3. The bipolaron binding energy W in cylindrical (a) and planar (b) quantum wires 
plotted versus the dimensionless radius R and width L, respectively. 

Fig. 4. The minimal value (at rj — 0) of the critical electron-phonon coupling constant a c 
plotted versus R and L in cylindrical (a) and planar (b) quantum wires, respectively. 

Fig. 5. The ratio of the critical Coulomb repulsion constant U c and a, as a function of a, in 
cylindrical quantum wires for R = 0.01 (1), 0.5 (2), 1.0 (3), and 20.0 (4). 

Fig. 6. The bipolaron effective mass {rn^ ip )^/m in cylindrical (a) and planar (b) quantum 
wires plotted versus R and L, respectively. The curves for the effective mass are broken off 
as the bipolaron becomes unstable. 



15 




Fig. 1 . A scheme of the trial system which contains two electrons connected with 
two "fictitious" particles through the elastic attraction and models the Coulomb 
interaction by the elastic repulsion. 




Fig. 2. A scheme of cylindrical and planar quantum wires. 
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Fig. 3. The bipolaron binding energy Win cylindrical (a) and 
planar (b) quantum wires plotted versus the dimensionless 
radius R and widthZ, respectively. 




Fig. 4. The minimal value (at T]=0) of the critical electron- 
phonon coupling constant oc c plotted versus R and L in 
cylindrical (a) and planar (b) quantum wires, respectively. 
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Fig. 5. The ratio of the critical Coulomb repulsion constant 
U c and oc, as a function of oc, in cylindrical quantum wires for 
£=0.01 (1), 0.5 (2), 1.0 (3), and20.0 (4). 




Fig. 6. The bipolaron effective mass m h Jm in cylindrical (a) 

and planar (b) quantum wires plotted versus R and L, 
respectively. The curves for the effective mass are broken off 
as the bipolaron becomes unstable. 



